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Abstract:  In  the  present  work,  a  ballistic  material-model  development  approach  is  combined 
with  transient  non-linear  dynamics  simulations  of  the  projectile/armour  interactions  to  explore 
the  armour-hard-facing  potential  of  multi-walled  carbon  nanotube  (MWCNT)  reinforced,  poly¬ 
vinyl-ester-epoxy  (PVEE) -matrix  composite  mats.  This  approach  is  applied  to  improving  the 
ballistic-protection  performance  of  E-glass  fibre-mat  reinforced  PVEE-matrix  laminate 
armour.  Two  different  architectures  of  the  MWCNT-reinforced/PVEE-matrix  composite  mats 
were  considered:  (a)  a  MWCNT-ply  mat  structure  in  which  the  MWCNT  reinforcements  are 
aligned  parallel  to  the  armour  faces  and  (b)  a  MWCNT-forest  mat  structure  in  which  the 
MWCNT  reinforcements  are  aligned  orthogonally  to  the  armour  faces.  The  projectile/armour 
interaction  simulation  results  showed  that,  at  low  volume  fractions  of  MWCNTs,  both 
armour  architectures  yield  no  discernable  increase  in  the  armour  ballistic-protection  perform¬ 
ance  as  measured  by  the  armour.  On  the  other  hand,  at  high  MWCNT  volume  fractions  of 
MWCNTs,  the  first  armour  architecture  remained  ineffective  whereas  the  second  showed  a 
minor  improvement  in  the  ballistic-protection  performance  relative  to  the  corresponding 
monolithic  armour.  These  results  were  rationalized  using  published  experimental  observations 
pertaining  to  the  effect  of  MWCNTs  on  the  in-plane  and  the  through- the-thickness  properties  of 
fibre-mat/polymer-matrix  composite  materials. 
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1  INTRODUCTION 

Recent  efforts  of  the  US  Army  have  been  aimed  at 
becoming  more  mobile,  deployable,  and  sustainable 
while  maintaining  or  surpassing  the  current  levels 
of  lethality  and  survivability.  Current  battlefield 
vehicles  have  reached  in  excess  of  70  tons  due  to 
ever  increasing  lethality  of  ballistic  threats,  which 
hinders  their  ability  to  be  readily  transported 
and  sustained.  Therefore,  a  number  of  research  and 
development  programs  are  under  way  to  engineer 
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light-weight,  highly  mobile,  transportable,  and 
lethal  battlefield  vehicles  with  a  target  weight  under 
20  tons.  To  attain  these  goals,  significant  advances 
are  needed  in  the  areas  of  light-weight  materials 
development  and  multi-functional  integration  of 
armour. 

The  earliest  reports  of  composite  light  armour 
developments  appear  in  the  late  1960’s,  when  Wilkins 
and  co-workers  [1-3]  demonstrated  the  potential 
of  a  hybrid  armour  consisting  of  a  hard  ceramic 
strike  face  and  a  thin  ductile  composite  backing.  The 
good  ballistic  performance  of  the  armour  was  attribu¬ 
ted  to  the  low-density,  high  hardness,  and  high  com¬ 
pressive  strength  of  the  ceramic  compound  with  a 
high  weight-normalized  strength  and  stiffness  com¬ 
posite  material.  The  first  successful  integration  of 
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this  armour  came  in  the  late  1980’s  where  a  glass- 
reinforced  polyester  hull  covered  with  hard-faced  cer¬ 
amic  titles  was  developed  to  replace  the  aluminum 
hull  of  a  Bradley  Infantry  Fighting  Vehicle  [4].  The 
resulting  vehicle  is  generally  referred  to  as  the  Com¬ 
posite  Infantry  Fighting  Vehicle  [5].  A  composite, 
ground-up  vehicle  development  program,  referred 
to  as  the  Composite  Armoured  Vehicle  (CAV)  pro¬ 
gram,  was  subsequently  implemented  [6],  To  meet 
the  stringent  weight  and  ballistic  performance 
requirements  of  the  CAV,  the  concept  of  multi¬ 
functional  polymer-matrix  composite-based  armour 
was  developed.  A  schematic  of  such  armour  is 
depicted  in  Fig.  1. 

Each  layer  in  the  multi-functional  polymer-matrix 
composite-based  armour  serves  a  specific  purpose, 
yet  combinations  of  layers  provide  role-sharing 
multi-functionality.  A  thin,  protective  polymer- 
matrix  composite  face  sheet  on  the  outside  of  the 
vehicle  serves  to  protect  the  ceramic-ballistic  tiles 
from  incidental  damage.  The  primary  role  of  the  cer¬ 
amic  tiles  is  to  break  up  and/or  erode  the  projectile 
upon  impact.  The  subsequent  rubber-based  layer  is 
utilized  to  improve  multi-hit  ballistic  performance 
of  the  armour.  The  following  thick-section  composite 
plate  serves  as  the  structural  support  for  the  vehicle,  a 
structural  backing  for  the  ballistic  tiles,  and  also  to 
assist  in  catching  of  the  projectile  remnants  and 
ceramic-tile  fragments  and  in  absorbing  their 
residual  kinetic  energy.  Finally,  a  fire-protective 
‘spall’  layer  of  phenolic  polymer  is  incorporated  on 
the  inner  surface  of  the  vehicle.  Other  layers  can  be 
incorporated  to  provide  additional  functionality, 
such  as  electro-magnetic  ground  planes,  signature 
control,  etc.  [4]. 

Although  the  development  of  the  multi-functional 
armour  played  an  important  role  in  understanding 
the  performance  of  composite  armour,  the 
areal-density  normalized  ballistic  protection  per¬ 
formance  was  found  to  fall  far  short  of  the  current 
requirements. 

In  the  ongoing  work  [7-10],  the  potential  of  next- 
generation  polymer-matrix  composites  reinforced 
with  carbon  nanotubes  as  advanced  light-weight 
armour  is  being  investigated  both  computationally 
and  experimentally.  It  was  initially  anticipated  that, 
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Fig.  1  A  schematic  of  the  multi-functional  composite 
armour 


due  to  a  high  hardness  of  the  nanotubes,  carbon- 
nanotube  reinforced  polymer-matrix  composite 
armour  may  be  very  effective  in  eroding/ fracturing 
the  projectile  while  a  high  strength  combined  with  a 
high  ductility  of  these  composite  materials  was 
expected  to  make  them  very  efficient  in  absorbing 
the  projectile’s  kinetic  energy. 

Before  the  ballistic  performance  of  such  armour 
could  be  assessed  computationally,  a  dynamic-material 
model  had  to  be  derived  for  the  nanotube  reinforced 
polymer-matrix  composite  material.  This  was  done  in 
a  series  of  papers  [7-10]  using  a  multi-length  scale 
approach.  Specifically:  (a)  molecular-level  computer 
simulations  were  used  to  determine  the  material  mech¬ 
anical  properties  of  ‘effective  fibres’  (multi-walled 
carbon  nanotubes  (MWCNTs)  surrounded  with  a 
nanometre  thick  layer  of  poly-vinyl-ester-epoxy 
(PVEE)  with  an  altered  confirmation)  [7];  (b)  micro¬ 
mechanics  type  calculations  were  employed  to  deter¬ 
mine  the  effective  mechanical  properties  of  MWCNT 
reinforced  PVEE  composite  mats  [8] ;  and  (c)  continuum 
level  homogenization  techniques  were  used  to  deter¬ 
mine  the  effective  mechanical  properties  of  E-glass 
reinforced  composites  with  MWCNT-doped  PVEE 
matrix  [9].  Once  the  required  material  models  were 
derived,  a  series  of  transient  non-linear  dynamics  simu¬ 
lations  of  the  impact  of  a  composite  armour  with  a  frag¬ 
ment  simulating  projectile  (FSP),  were  carried  out  [10] 
(It  should  be  noted  that  the  term  FSP  is  used  to 
denote  a  simulated  projectile,  which  mimics  a  typical 
fragment  generated  during  explosion  of  an  improvised 
explosive  device  (IED).  In  all  of  the  cases  analysed  by 
the  present  authors,  in  this  work  and  in  prior  work,  ero¬ 
sion,  but  not  fragmentation  of  the  FSP  was  observed) . 

Two  distinct  architectures  of  the  armour  were  con¬ 
sidered:  (a)  hybrid  armour  consisting  of  a  100-|xm 
thick  MWCNT  reinforced  PVEE-matrix  mat  sand¬ 
wiched  between  two  E-glass  mat  reinforced  PVEE- 
matrix  lamina;  and  (b)  a  monolithic  E-glass  reinforced 
armour  in  which  the  PVEE  matrix  was  doped  with  the 
same  total  volume  of  MWCNTs  as  in  the  case  of  the 
hybrid  armour.  The  results  obtained  showed  that  the 
use  of  either  of  these  options  gives  rise  to  a 
relatively  small  (~6  per  cent)  increase  in  the  ballistic- 
protection  performance  (as  quantified  by  V50,  the 
initial  velocity  and  a  given  FSP  that  yields  a  50  per 
cent  probability  of  penetrating  the  armour)  of  the 
composite  armour.  This  finding  was  rationalized  by 
the  fact  that,  due  to  the  nature  of  the  composite- 
laminate  manufacturing  method  used,  MWCNT 
reinforcements  were  aligned  parallel  to  the  armour 
faces.  This  caused  only  a  minor  improvement 
(less  than  a  per  cent)  in  the  through-the-thickness 
properties. 

It  was  further  suggested  that  the  use  of  MWCNT 
reinforcements  aligned  orthogonally  to  the  armour 
faces  (‘MWCNT  forest’)  may  be  more  beneficial  to 
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the  improvement  of  the  armour’s  ballistic  protection 
performance.  This  recommendation  was  supported 
by  the  recent  work  of  Veedu  et  al.  [11],  reporting 
substantial  (in  excess  of  100  per  cent)  increase  of 
composite  laminate  properties  (in  particular,  the 
through-the-thickness  mechanical  properties)  in 
which  MWCNTs  were  aligned  in  the  direction 
normal  to  the  faces  of  the  composite-laminate. 

The  objective  of  the  present  work  is  two-fold:  (a)  to 
derive  and  parameterize  a  material  model  for 
MWCNT  forest  reinforced  PVEE  matrix  mats  and  (b) 
to  assess  the  potential  of  such  mats  used  as  armour 
hard-facing  to  improve  the  ballistic  protection  per¬ 
formance  of  E-glass  mat  reinforced  PVEE  matrix 
armour.  The  objective  (b)  will  be  addressed  in  the 
present  work  using  a  series  of  transient  non-linear 
dynamics  computational  analyses.  In  the  ongoing 
work,  the  objective  (b)  is  also  being  addressed  exper¬ 
imentally.  Part  of  the  ongoing  work  includes  fabrica¬ 
tion  of  the  MWCNT  forest-mat  hard-face.  These 
MWCNTs  aligned  in  a  perpendicular  direction  with 
respect  to  the  substrate  (a  plain-woven  mat  of 
E-glass  in  the  present  work)  are  grown  in  a  tubular 
furnace  using  a  chemical  vapour  deposition  (CVD) 
process.  The  plain-woven  E-glass  fabric  is  used  as 
an  internal  lining  of  the  furnace.  A  reaction-gas  mix¬ 
ture  (consisting  of  a  carbon- containing  pre-cursor 
gas,  a  catalytic -particles  pre-cursor  gas,  and  a  carrier 
gas)  is  fed  into  the  furnace  whose  (axial)  temperature 
profile  is  adjusted  to  yield  a  uniform  deposition  rate 
of  MWCNTs  along  the  length  of  the  furnace.  The 
MWCNTs  grow  orthogonally  to  the  local  surface  of 
the  furnace  inner  lining  whereas  the  catalytic  par¬ 
ticles  reside  on  the  top  (growing)  end  of  the 
MWCNTs.  This  process  can  yield  MWCNT  forests 
with  an  average  MWCNT  length  of  150-250  pm. 
Once  the  MWCNT -forest  coated  E- glass  mat  furnace 
lining  is  removed  and  placed  on  top  of  the  E-glass 
mats  (with  the  MWCNT-forest  side  facing  outward), 
a  standard  vacuum  assisted  resin  transfer  moulding 
(VARTM)  process  is  used  to  produce  the  composite 
armour  with  a  MWCNT-forest  hard-face.  Upon  the 
completion  of  fabrication  of  the  hard-faced  armour, 
its  ballistic  performance  is  tested  by  striking  armour 
with  cylindrical  steel  projectiles  of  differing  masses 
and  impact  velocities. 

The  organization  of  the  paper  is  as  follows.  A  brief 
overview  of  the  non-linear  dynamics  computational 
procedure  utilized  in  the  present  work  is  given  in  sec¬ 
tion  2.1.  Constitutive  models  used  to  represent  the 
behaviour  of  the  projectile  and  the  armour  materials 
under  ballistic  impact  conditions  are  discussed  in 
section  2.2.  Details  of  the  numerical  model  used  to 
analyse  the  impact  and  penetration  of  the  armour 
by  a  fragment  simulating  projectile  are  presented  in 
section  2.3.  The  results  obtained  in  the  present  work 
are  presented  and  discussed  in  section  3.  The  main 


conclusions  resulting  from  the  present  work  are  sum¬ 
marized  in  section  4. 


2  COMPUTATIONAL  PROCEDURE 

2.1  Transient  non-linear  dynamics  modelling  of 
projectile/target  interactions 

Computational  assessment  of  the  ballistic-protection 
performance  of  armour  with  respect  to  IED  frag¬ 
ments  is  generally  done  by  carrying  out  numerical 
simulations  of  the  interactions  between  an  FSP  and 
armour  (an  E- glass  mat  reinforced  PVEE  matrix  com¬ 
posite  laminate  hard-faced  with  a  MWCNT-forest- 
based  mat,  in  the  present  work).  Such  calculations 
have  been  performed,  in  the  present  work,  using 
AUTODYN,  a  general-purpose  transient  non-linear 
dynamics  modelling  and  simulation  software  [12]. 
Within  AUTODYN,  the  appropriate  mass,  momen¬ 
tum  and  energy  conservation  equations  are  com¬ 
bined  with  the  attendant-materials  model  equations 
and  the  appropriate  initial  and  boundary  conditions 
and  solved  numerically  using  a  second- order  accu¬ 
rate  explicit  scheme.  The  numerical  framework  (i.e. 
the  ‘processor’  as  referred  to  in  AUTODYN)  used  is 
generally  dependent  on  the  physical  nature  of  the 
problem  being  studied  and,  for  multi-domain  pro¬ 
blems,  different  domains  can  be  analysed  using 
different  processors.  The  Lagrange  processor  or  the 
SPH  (smooth  particle  hydrodynamics)  processor  are 
typically  used  for  solid-continuum  structures, 
whereas  the  Euler  processor  is  commonly  used  for 
modelling  gases,  liquids,  or  solids  undergoing  large 
deformations  and  density  changes.  The  Lagrange- 
based  shell  and  beam  processors  are  designated  for 
modelling  shell-  and  beam-like  solid  structures, 
respectively. 

In  the  present  work,  the  ballistic-protection  per¬ 
formance  of  E-glass  mat  reinforced  polymer-matrix 
composite  armours  hard-faced  with  a  MWCNT - 
forest-based  mat  under  FSP  threats  was  analysed 
using  the  Lagrange  processor.  The  interactions 
between  the  projectile  and  the  armour  are  accounted 
for  through  the  use  of  the  sub-domain  interaction 
options  within  AUTODYN  [12],  which  were  over¬ 
viewed  in  detail  in  our  recent  work  [10].  Also,  a 
detailed  discussion  regarding  the  effect  of  the  pro¬ 
cessor  choice  (Lagrange  versus  SPH)  for  the  projectile 
and  the  target  on  the  computational  results  can  be 
found  in  reference  [8] . 


2.2  Material  models 

The  complete  formulation  of  a  transient  non-linear 
dynamics  problem  involving  impact  of  armour  by 


JMDA136  ©  IMechE  2008 


Proc.  IMechE  Vol.  222  Part  L:  J.  Materials:  Design  and  Applications 


18 


M  Grujiclc,  W  C  Bell,  S  B  Biggers,  K  L  Koudela,  J  F  Tarter,  and  B  A  Cheeseman 


an  FSP  entails  the  knowledge  of  materials  models 
(material-specific  relations  between  pressure,  stress, 
mass  density,  strain,  strain  rate,  internal  energy  den¬ 
sity,  etc.).  These  relations  typically  involve:  (a)  an 
equation  of  state;  (b)  a  strength  equation;  (c)  a  failure 
equation;  and  (d)  an  erosion  equation  for  each  con¬ 
stituent  material.  The  equation  of  state  defines 
pressure  dependence  on  mass  density  and 
internal-energy  density  (and,  in  the  case  of  aniso¬ 
tropic  materials,  on  deviatoric  strain).  The  strength 
and  failure  equations  define  the  evolutions  of  the 
deviatoric  stress  in  the  elastic  regime,  elastic -plastic 
regime,  and  in  the  post  failure  initiation  regime.  In 
other  words,  the  equation  of  state  along  with  the 
strength  and  failure  equations  enable  assessment 
of  the  evolution  of  the  complete  stress  tensor 
during  a  transient  non-linear  dynamics  analysis. 
Such  an  assessment  is  needed  where  the  governing 
conservation  equations  are  being  solved.  The  erosion 
equation  is  generally  intended  for  eliminating 
numerical  solution  difficulties  arising  from  highly 
disordered  Lagrangian  cells.  Nevertheless,  the  ero¬ 
sion  equation  is  often  used  to  provide  additional 
material  failure  mechanism  especially  in  materials 
with  limited  ductility.  When  a  computational/ 
material  cell  is  eroded,  the  freed  nodes  are  retained 
along  with  their  velocities  in  order  to  conserve 
momentum  of  the  system. 

In  the  following  two  sections,  a  brief  description  is 
given  of  two  material  models  used  in  the  present 
work.  The  first  (isotropic)  material  model  is  used  to 
represent  the  behaviour  of  an  A1S1  4340-steel  FSP. 
The  values  of  the  material  parameters  for  AISI  4340 
steel  are  taken  from  the  AUTODYN  materials  library 
[12].  The  other  (transversely  isotropic)  material 
model  is  used  to  represent  the  behaviour  of  the 
E-glass  mat  reinforced  MWCNT-doped  PVEE  compo¬ 
site  laminates  as  well  as  the  MWCNT-reinforced  PVEE 
composite  mats  (for  both  the  case  of  MWCNTs 
aligned  parallel  with  and  orthogonal  to  the  mat 
faces).  The  corresponding  material  parameters  for 
the  E-glass  mat  reinforced  MWCNT-doped  PVEE 
composite  laminates  and  the  MWCNT-reinforced 
PVEE  composite  mats  with  the  MWCNTs  aligned 
orthogonal  to  the  mat  faces  are  taken  from  the  pre¬ 
vious  work  [10].  A  derivation  of  the  material  par¬ 
ameters  for  the  composite  material  with  the 
MWCNTs  aligned  normal  to  the  mat  faces  is  carried 
out  in  section  3.1  in  the  present  work. 

2.2.1  AISI  4340  steel 

The  material  model  for  AISI  4340  steel  includes:  (a)  an 
(internal -energy  density  invariant)  linear  equation  of 
state;  (b)  a  strength  modes  consisting  of  the  von 
Mises  yield  criterion,  the  Prandtl-Reuss  associated 
flow  rule,  and  the  Johnson-Cook  flow-strength 


model;  (c)  a  stress-based  damage-initiation  isotropic¬ 
softening  ductile-failure  Johnson-Cook  model;  and 
(d)  an  erosion  model  based  on  an  instantaneous  geo¬ 
metrical  strain  of  2.0.  Since  a  detailed  account  of  the 
material  model  for  AISI  4340  steel  was  given  in  the 
recent  work  [9] ,  no  further  details  of  this  model  will 
be  presented  here. 

2.2.2  Planar  isotropic  material  model  for  E-glass 
mat  and  MWCNT  reinforced  poly-vinyl-ester 
epoxy  composite  laminates 

As  discussed  earlier,  the  same  planar  isotropic 
material  model  (presented  in  this  section)  is  used 
(with  different  parameterizations)  to  describe  the 
behaviour  of  five  materials: 

(a)  a  composite  material  consisting  of  a  0.49-vol% 
MWCNT-doped  PVEE  matrix,  which  is  reinforced 
with  53  vol%  E-glass  mats; 

(b)  a  composite  material  consisting  of  a  0.49-vol% 
MWCNT-doped  PVEE  matrix,  which  is  reinforced 
with  10  vol%  MWCNT  mats  aligned  parallel  with 
the  armour  strike  face; 

(c)  a  composite  material  consisting  of  a  0.49-vol% 
MWCNT-doped  PVEE  matrix  which  is  reinforced 
with  30  vol%  MWCNT  mats  aligned  parallel  with 
the  armour  strike  face; 

(d)  a  composite  material  consisting  of  a  0.49-vol% 
MWCNT-doped  PVEE  matrix,  which  is  reinforced 
with  10  vol%  MWCNT-forest  mats; 

(e)  a  composite  material  consisting  of  a  0.49-vol% 
MWCNT-doped  PVEE  matrix,  which  is  reinforced 
with  30  vol%  MWCNT-forest  mats. 

It  should  be  noted  that  material  (a)  is  used  to  con¬ 
struct  the  main  body  of  armour  whereas  materials 
(b)  -  (e)  are  used  (one  at  a  time)  as  100-|jim  thick  hard- 
face.  For  brevity,  the  model  used  to  describe  mechan¬ 
ical  behaviour  of  the  five  materials  given  above  (with 
different  parameterizations)  will  be  referred  to  as  the 
composite-laminate  model  in  the  remainder  of  the 
manuscript.  Before  this  material  model  is  presented, 
it  should  be  noted  that  the  materials  described  are  all 
planar-isotropic  with  an  x2~x3  plane  of  isotropy  and 
the  through-the-thickness  direction  parallel  with  the 
Xi-axis.  Also,  since  a  detailed  account  of  this  material 
model  was  given  in  the  recent  work  [10],  only  a  brief 
overview  of  the  key  relations  will  be  given  in  this 
section. 

2.2.2. 1  Equation  of  state.  A  polynomial  form  of 
the  equation  of  state  is  used  to  represent  the  pressure 
dependence  on  mass  density  and  the  internal 
energy  density  (i.e.  temperature).  Furthermore, 
due  to  the  non-isotropic  nature  of  the  material  at 
hand,  an  additional  term  is  used  to  represent  the 
pressure  dependence  of  the  deviatoric  strain 
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components.  Such  a  term  is  not  present  in  isotropic 
materials. 

2.2.2.2  Strength  model.  The  evolution  of  the  devia- 
toric  part  of  the  stress  tensor  is  described  using  the 
planar-isotropic  elastic  stiffness  tensor,  a  six- 
parameter  parabolic  yield  function,  the  Prandtl- 
Reuss  normality  flow  rule  and  a  piece-wise  linear 
strain  hardening  constitutive  material  flow-strength 
equation.  The  elastic  stiffness  tensor  is  used  for 
deviatoric-stress  calculations  during  elastic  loading/ 
unloading  and  elastic  reloading,  whereas  the  yield 
criterion,  the  flow  rule,  and  the  flow-strength  relation 
are  used  to  update  the  stress  during  and  elastic-plas¬ 
tic  deformation  step. 

2.2.2.3  Failure  model.  Once  the  material  reaches 
the  condition  for  damage  initiation,  the  stress  evol¬ 
ution  is  controlled  by  the  progression  of  damage 
and  not  by  plastic  deformation.  The  failure  model 
used  here  utilizes  a  quadratic  stress-based  failure 
initiation  criterion  (one  equation  for  each  of  the 
three  principal  material  planes).  Consequently,  the 
failure  initiation  entails  the  knowledge  of  six  (three 
normal  and  three  shear)  failure-initiation  stresses. 
Once  the  failure  initiation  has  occurred  on  a  given 
principal  plane,  the  material  begins  to  accumulate 
damage  which  in  turn,  causes  the  quadratic  failure 
surface  to  shrink.  The  resulting  damage-induced  irre¬ 
versible  deformation  is  also  governed  by  the  Prandtl- 
Reuss  normality  flow  rule.  The  ultimate  failure  on  a 
give  principle  plane  takes  place  when  the  combined 
damage  arising  from  the  corresponding  normal  and 
two  shear  stress  components  acquire  its  maximum 
allowed  value.  Such  a  value  is  controlled  by  the  corre¬ 
sponding  fracture  energies.  Thus,  the  damage  evol¬ 
ution  part  of  the  failure  model  entails  the 
knowledge  of  six  (three  normal  and  three  shear)  frac¬ 
ture  energies. 

2.2.2A  Erosion  model.  As  discussed  in  the  pre¬ 
vious  work  [8],  the  computational  results  attained 
during  a  transient  non-linear  dynamics  simulation 
of  FSP/ armour  interactions  are  quite  sensitive  to 
the  choice  of  erosion  equation  parameters,  where 
the  FSP  and  the  armour  are  both  represented  using 
Lagrange  models.  Consequently,  a  substantial  effort 
was  invested  in  the  previous  work  in  order  to  arrive 
at  the  value  of  the  erosion  strain,  which  ensures 
that  both  the  ballistic  protection  (as  quantified  by 
V50),  and  the  distribution  of  the  damage  as  predicted 
by  the  model  are  in  good  agreement  with  the  exper¬ 
imental  counterparts.  A  strain  value  of  0.65  was 
found  to  be  the  optimal  value  for  material  (a)  and 
0.70  for  materials  (b)-(e). 


2.3  Problem  definition  and  computational 
analysis 

As  stated  earlier,  one  of  the  objectives  of  the  present 
work  is  to  carry  out  transient  non-linear  dynamics 
analyses  of  the  impact  and  penetration  of  E-glass 
mat  reinforced  PVEE  matrix  composite  armour  hard- 
faced  with  a  MWCNT-reinforced  PVEE-matrix  mat  by 
an  FSP  in  order  to  determine  the  ballistic-protection 
performance  of  such  armour.  Three  distinct  armour 
architectures  were  considered:  (a)  a  12.9-mm  thick 
monolithic  composite  armour  based  on  E-glass 
reinforced  mats  and  a  0.49-vol%  MWCNT-doped 
PVEE  matrix;  (b)  a  hybrid  armour  based  on  a 
0.2-mm  thick  10  or  30vol%  MWCNT-forest 
reinforced  PVEE-matrix  hard-face  mat  and  a 
12.7-mm  thick  backing  made  of  the  same  E-glass 
reinforced  mats  and  a  0.49-vol%  MWCNT-doped 
PVEE  matrix,  as  in  (a);  and  (c)  the  same  as  in  (b) 
but  with  the  MWCNT  reinforcements  parallel  with 
the  armour  faces.  A  summary  of  the  five  armour 
architectures  is  provided  in  Tables  1  to  3. 

The  work  was  limited  to  the  case  of  a  normal 
impact  of  the  armour  by  the  FSP  and,  due  to  the 
attendant  axi-symmetric  nature  of  the  problem,  all 
the  calculations  were  carried  out  using  a  two- 
dimensional  (axisymmetric)  model.  A  simple  sche¬ 
matic  of  the  projectile /armour  impact/penetration 
problem  analysed  here  is  given  in  Fig.  2.  The  (0.30 
caliber  FSP)  projectile  is  cylindrical  in  shape  with  a 
diameter  of  7.62  mm  and  a  height  of  7.62  mm.  The 
initial  normal  velocity  of  FSP  was  varied  in  a  range 
between  400  and  500  m/s.  A  fixed  overall  thickness 
of  12.9  mm  and  a  fixed  lateral  dimension  of  80  mm 
for  the  armour  are  used  for  all  the  cases  analysed. 

The  Lagrange  processor  was  used  to  represent  both 
the  FSP  and  the  armour.  The  projectile  was  analysed 
using  a  mesh  of  200  rectangular  cells  and  an  armour 
consisting  of  540  cells.  To  improve  the  accuracy  of 
the  analysis,  smaller  cells  were  used  in  the  regions  of 
the  projectile  and  the  armour  involved  in  the  projec¬ 
tile/armour  interactions.  The  interaction  between  the 
FSP  and  the  armour  was  accounted  for  using  the 
parts  coupling  option  available  in  AUTODYN  [12].  A 
constant  friction  coefficient  of  0.1  between  FSP  and 
armour  is  used  through  all  the  calculations  carried 
out  in  the  present  work.  Except  for  the  projectile/ 
armour  contact  surfaces  and  lateral  faces  of  the 
armour,  (default)  zero-stress  boundary  conditions 
were  prescribed  on  all  faces  of  the  projectile  and  the 
armour.  To  ensure  that  the  blast  waves  are  not  reflected 
at  the  lateral  faces  of  the  armour,  the  transmit-type 
boundary  conditions  were  applied  on  the  faces.  A  stan¬ 
dard  mesh  sensitivity  analysis  was  conducted  (the 
results  not  shown  for  brevity)  in  order  to  ensure  that 
the  results  obtained  were  effectively  insensitive  to  vari¬ 
ations  in  the  size  of  the  cells  used. 
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Table  1  Orthotropic  equation  of  state  model  parameters  for  materials  used  in  the  present  work 


Value 


MWCNT  doped  derakane  510A  MWCNT  doped  derakane 
plus  in-plane  MWCNT  mat  510A  plus  MWCNT  forest  mat 
MWCNT  doped  E-glass  -  - 


Parameter 

Symbol 

Unit 

derakane  510A 

10% 

30% 

10% 

30% 

Density 

P 

g/cm3 

1.9349 

1.23 

1.23 

1.23 

1.23 

Young’s  modulus  1 1 

Eii 

kPa 

2.0899 

X 

107 

3.2320  x  107 

6.1126 

x  106 

6.1126  x  106 

8.7189  x  107 

Young’s  modulus  22 

E22 

kPa 

4.3046 

X 

107 

4.6053  x  106 

1.3806 

x  107 

1.3806  x  107 

5.9408  x  106 

Young’s  modulus  33 

E33 

kPa 

4.3046 

X 

107 

4.6053  x  106 

1.3806 

x  107 

1.3806  x  107 

5.9408  x  106 

Poisson’s  ratio  12 

V12 

- 

0.1751 

0.353  714 

0.174  786 

0.174  786 

0.293  245 

Poisson’s  ratio  23 

^23 

- 

0.1080 

0.598  707 

0.342  873 

0.342  873 

0.293  245 

Poisson’s  ratio  31 

14 1 

- 

0.3608 

0.050  402 

0.394  793 

0.394  793 

0.019  981 

Shear  modulus  12 

G12 

kPa 

6.0290 

X 

10s 

1.3776  x  106 

1.4086 

x  10s 

1.4086  x  106 

1.7500  x  106 

Shear  modulus  23 

G23 

kPa 

6.7433 

X 

10s 

1.4404  x  106 

3.1847 

x  10s 

3.1847  x  10s 

2.2000  x  106 

Shear  modulus  31 

G3I 

kPa 

6.0290 

X 

106 

1.3776  x  106 

1.4086 

x  106 

1.4086  x  106 

1.7500  x  106 

Bulk  modulus 

K 

kPa 

1.9579 

X 

107 

8.7204  x  106 

8.8361 

x  107 

8.8361  x  107 

1.5762  x  107 

Volumetric  parameter 

^2 

kPa 

5.0000 

X 

107 

5.0000  x  107 

5.0000 

x  107 

5.0000  x  107 

5.0000  x  107 

Volumetric  parameter 

^3 

kPa 

0 

0 

0 

0 

0 

Volumetric  parameter 

B0 

kPa 

0 

0 

0 

0 

0 

Volumetric  parameter 

Bi 

kPa 

0 

0 

0 

0 

0 

Specific  heat 

Cy 

J/kgK 

1.4200 

X 

103 

1.4200  x  103 

1.4200 

x  103 

1.4200  x  103 

1.4200  x  103 

Reference  temperature 

Tref 

K 

300 

300 

300 

300 

300 

Table  2  Orthotropic  yield  strength  model  parameters  for  materials  used  in  the  present  work 


Value 


Parameter 

Symbol 

Unit 

MWCNT  doped  E-glass 
derakane  510A 

MWCNT  doped  derakane  510A  plus 
in-plane  MWCNT  mat 

MWCNT  doped  derakane  510A  plus 
MWCNT  forest  mat 

10% 

30% 

10% 

30% 

Stress  parameter  1 1 

^11 

_ 

40.5730 

0.1038 

1.1107 

1.1107 

0.0200 

Stress  parameter  22 

A2  2 

- 

1.0000 

1.0000 

1.0000 

1.0000 

1.0000 

Stress  parameter  33 

A33 

- 

1.0000 

1.0000 

1.0000 

1.0000 

1.0000 

Stress  parameter  12 

a12 

- 

-8.1789 

-0.3646 

-0.5395 

-0.5395 

-0.4522 

Stress  parameter  23 

A23 

- 

-8.1789 

-0.3646 

-0.5395 

-0.5395 

-0.4522 

Stress  parameter  31 

7131 

- 

-0.5000 

-0.5000 

-0.5000 

-0.5000 

-0.500 

Stress  parameter  44 

G12 

- 

1.5000 

1.500 

1.5000 

1.5000 

1.500 

Stress  parameter  55 

G23 

- 

28.9650 

0.916  570 

1.5950 

1.5950 

0.9622 

Stress  parameter  66 

G31 

- 

28.9650 

0.916  570 

1.5950 

1.5950 

0.9622 

Effective  stress 

°Eff 

kPa 

7.7434  x  105 

6.4924  x  104 

8.8678  x  104 

8.8678  x  104 

1.6754  x  106 

Table  3  Orthotropic  softening  failure  model  parameters  for  materials  used  in  the  present  work 


Value 


MWCNT  doped  derakane  510A  MWCNT  doped  derakane  510A 
plus  in-plane  MWCNT  mat  plus  MWCNT  forest  mat 

MWCNT  doped  E-glass  -  - 


Parameter 

Symbol 

Unit 

derakane  510 

10% 

30% 

10% 

30% 

Failure  stress 

on 

kPa 

1.19110  x  105 

8.2441 

X 

104 

9.7734 

X 

104 

1.7271 

X 

105 

8.0893  x  105 

Failure  stress 

°22 

kPa 

6.5121  x  105 

8.6886 

X 

104 

1.1674 

X 

105 

6.3612 

X 

104 

2.2618  x  105 

Failure  stress 

U33 

kPa 

6.5121  x  105 

8.6886 

X 

104 

1.1674 

X 

105 

6.3612 

X 

104 

2.2618  x  105 

Failure  stress 

<Ti2 

kPa 

9.9680  x  104 

4.8648 

X 

104 

6.0431 

X 

104 

4.6983 

X 

104 

1.8770  x  105 

Failure  stress 

U23 

kPa 

3.7598  x  105 

5.0164 

X 

104 

6.7400 

X 

104 

3.6727 

X 

104 

1.6821  x  105 

Failure  stress 

ff31 

kPa 

9.9680  x  104 

4.8648 

X 

104 

6.0431 

X 

104 

4.6983 

X 

104 

1.8770  x  105 

Fracture  energy 

G„ 

J/m2 

744 

927 

1100 

864 

2045 

Fracture  energy 

G22 

J/m2 

326 

1086 

1459 

716 

745 

Fracture  energy 

G33 

J/m2 

326 

1086 

1459 

716 

745 

Fracture  energy 

G12 

J/m2 

1246 

730 

906 

705 

716 

Fracture  energy 

G23 

J/m2 

846 

752 

1011 

551 

573 

Fracture  energy 

G31 

J/m2 

1246 

730 

906 

705 

716 

Erosion  strain 

e 

- 

0.65 

0.70 

0.70 

0.70 

0.70 
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Fig.  2  An  example  of  the  computational  domain  used 
in  the  present  work  to  analyse  the  interactions 
between  a  FSP  and  composite  armour  hard- 
faced  with  a  MWCNT-forest  mat 

3  RESULTS  AND  DISCUSSION 

3.1  Material  model  for  MWCNT-forest  reinforced 
PVEE  matrix  mats 

As  discussed  earlier,  the  material  model  for  E-glass 
mat-reinforced  MWCNT-doped  PVEE-matrix  com¬ 
posite  and  for  the  MWCNT-reinforced  PVEE-matrix 
mats  (in  which  the  MWCNT  reinforced  mats  are 
aligned  parallel  with  the  mat  faces)  were  determined 
in  reference  [10].  In  this  section,  a  brief  description  is 
given  of  the  procedure  used  to  determine  the  material 
model  for  MWCNT-forest  reinforced  PVEE-matrix 
mats.  Before  this  procedure  is  presented,  however,  a 
brief  description  is  given  of  the  process,  which  is 
being  utilized  in  the  ongoing  work  to  process  such 
mats. 


3.2  Derivation  of  material  model  for 

MWCNT-forest  reinforced  PVEE-matrix 
composite  mats 

In  this  section,  an  attempt  is  made  to  derive  a  simple 
material  model  for  the  MWCNT-forest  reinforced 
PVEE-matrix  composite  mats.  The  model  developed 
in  the  section  is  based  on  three  assumptions:  (a) 
the  in-plane  composite-material  properties  are 
controlled  by  the  PVEE-matrix;  (b)  the  through-the- 
thickness  composite  material  properties  are  con¬ 
trolled  by  the  MWCNT-forest;  and  (c)  the  effect  of 
the  interactions  between  the  PVEE  matrix  with  the 
MWCNT  forest  can  be  accounted  for  using  the 
recently  developed  planar-isotropic  material  model 
for  the  MWCNT-reinforced  PVEE-matrix  composite 
[7-8].  A  summary  of  the  material  model  parameters 
for  the  equation  of  state,  strength  model,  failure 
model,  and  the  erosion  model  for  (10  and  30  vol%) 
MWCNT-forest  reinforced  PVEE-matrix  composite 
mats  obtained  in  the  present  work  are  given  in 
Tables  1  to  3,  respectively.  It  should  be  noted  that 
within  the  coordinate  system  chosen  in  the  present 


work,  the  through-the-mat-thickness  direction 
coincides  with  the  Xi-axis. 

It  should  be  observed  that  Tables  1  to  3  contain  a 
large  number  of  parameters.  These  parameters  were 
obtained  in  the  previous  work  [7,  8]  by  conducting 
molecular-level  mechanical-test  analyses  with 
micromechanics  homogenization  procedures.  The 
results  of  these  procedures  were  then  fitted  using 
the  composite  laminate  material  model  [12],  which 
contains  the  equations  of  state  parameters  (A2,  A3, 
B0,  Bi),  strength  parameters  (An  to  A31,  Gi2,  G23, 
G31),  and  failure  parameters  (o-u  to  <r31,  Gn  to  G13). 
A  complete  experimental  validation  of  the  material 
model  was  not  carried  out,  but  the  results  for  V50 
and  the  extent  of  surrounding  damage  observed  in 
ballistic  tests  were,  in  the  case  of  in-plane 
MWCNT-mats  reinforced  PVEE-matrix  composite, 
predicted  be  the  model. 

It  should  also  be  noted  that  no  direct  experimental 
validation  of  the  material  model  was  carried  out  in 
the  present  work.  The  sponsor  of  this  project  was 
interested  in  a  physics-based  model  (obtained  by 
modifying  an  existing  model  for  polymer-matrix 
laminates),  which  can  provide  reasonable  predictions 
of  V50  and  the  extent  of  surrounding  damage  in  the 
armour.  The  ongoing  experimental  work  [13]  has 
yielded  preliminary  results,  which  suggest  that  the 
model  developed  here  provides  a  reasonable  predic¬ 
tion  of  armour  ballistic-impact  conditions. 

3.2.1  Equation  of  state 

Due  to  the  planar-isotropic  nature  of  the  material  at 
hand,  an  orthotropic  material  model  (with  interchan¬ 
geability  of  the  x2  and  x3  direction  related  properties) 
is  developed.  Within  this  model  the  pressure,  P, 
dependence  on  the  compaction,  p=(p/p0-  1), 
(where  p  is  density  and  p0  initial  density)  and  temp¬ 
erature,  T,  is  defined  as 

P  —  -K/jl  +  A2ijl2  -  A3p3  +  (B0  -  Bin)p0e 

—  g  (Cii  +  2C|2)ef| 

—  2  (C12  +  C22  +  C23)(e22  +  £33)  (1) 

where  K  is  the  bulk  modulus  defined  as 

K  —  -[Cn  +  2C22  +  4Ci2  +  2C23]  (2) 

where  e=  CV(T  —  Tn.f)  is  the  internal  energy  density, 
Cv  is  the  material  specific  heat,  Tre{  is  an  inertial 
energy- density  reference  temperature,  Q/s  are  the 
material  stiffness  coefficients  (coefficients  of  the 
material  6x6  stiffness  matrix),  effs  are  the  com¬ 
ponents  of  the  deviatoric  strain  matrix,  and  A2,  A3, 
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B0,  Bi  are  material-specific  parameters.  The  last  term 
on  the  right-hand  side  of  equation  (1)  represents  the 
coupling  between  pressure  and  the  deviatoric  strain 
and  is  absent  in  isotropic  materials. 

Equation  (1)  shows  that  a  complete  definition  of 
the  equation  of  state  entails  the  knowledge  of  the 
planar-isotropic  elastic-constants  of  the  material 
(Cn,  C-z2,  C12,  C23,  C44,  and  C55)  as  well  as  the  specifi¬ 
cation  ofA2,  A3,  B0,  Bi}  Cy,  and  Tre f.  It  should  be  noted 
that  only  normal  deviatoric  strain  components 
appear  in  the  pressure-deviatoric  strain  coupling 
term.  This  is  not  an  approximation  but  a  direct 
result  of  a  simple  derivation  of  the  pressure- 
deviatoric  strain  coupling  [12]. 

The  aforementioned  elements  of  the  elastic  stiffness 
matrix  were  determined  using  the  micromechanics - 
based  Mori-Tanka  type  elastic  homogeneity  model 
developed  in  the  previous  work  [8] . 

Due  to  the  similarity  in  the  chemical  make-up  and 
microstructure  between  the  MWCNT-forest  mats  and 
the  ones  in  which  the  MWCNT  reinforcements  are 
parallel  with  the  mat  faces  (the  model  for  the  latter 
material  was  derived  in  the  previous  work  [10]),  the 
values  for  A2,  A3,  B0,  Bt  Cv,  and  Tref  are  set  equal  in 
the  two  materials. 


3.2.2  Strength  model 

The  strength  model  is  defined  as  a  combination  of 
the  yield  criterion,  the  plastic -flow  rule  and  the 
constitutive  material  model.  To  account  for  the 
planar-isotropic  nature  of  the  material  at  hand,  a  six- 
parameter  parabolic  yield-function  (criterion)  is  used 
which  can  be  written  as 


d\\ 


°22,y 

°ll,y 


(5) 


^12 


k 


°"l2,ybal_biaxial 


(6) 


a2  3  =  —  0.5a22 


(7) 


(8) 


a55  —  ial  1  —  a22  ~  2Ul2) 


O) 


It  should  be  noted  that  some  of  the  relations 
between  different  a^’s  (e.g.  equation  (6))  were 
obtained  from  the  recognition  of  the  material  iso¬ 
tropy  in  the  x2-x3  plane.  Consequently,  the  corre¬ 
sponding  relations  in  the  isotropic  materials  were 
used. 

The  direction  of  the  plastic  flow  is  defined  using  the 
Prandtl-Reuss  normality  flow  rule.  As  far  as  the  con¬ 
stitutive  equation  (the  equation,  which  defines  the 
rate  at  which  material  flow-strength  increases  as  the 
equivalent  plastic  strain  increases)  is  concerned,  a 
piece-wise  linear  form  is  used.  It  should  be  noted 
that  an  analogous  procedure  for  the  derivation  of  a 
material  model  for  the  PVEE  matrix  composite 
material  reinforced  with  MWCNTs  aligned  parallel 
with  the  mat  faces  was  used  in  the  recent  work  [10], 
where  additional  details  for  the  procedure  can  be 
found. 


—  Anon  +  ^22(^22  T  ^33) 

+  2ai2crii(u22  +  033)  +  2a23(T22  0-33 
2^44(723  +  2r?55  T  cr^2 )  —  k  <  0  (3) 

Among  the  seven  material- specific  parameters 
^11-^66  and  k,  one  can  be  set  independently.  In  the 
present  work,  a22  is  set  to  1.  It  should  be  noted  that 
only  nine  out  of  ten  parameters  appearing  in 
equation  (3)  are  independent,  hence  one  of  these  par¬ 
ameters  can  be  set  independently  without  any  effect 
on  the  computational  results. 

The  values  of  the  remaining  parameters  in 
equation  (5)  are  defined  using  the  yield  stress  values 
corresponding  to  the  x1(  x2  (o-lliy  and  o-22,y,  respect¬ 
ively,)  and  balanced-biaxial  x1-x2  yield  stress 

1  U"l  2  ,y,bal_biaxial) 
k  —  °22,y 


3.2.3  Failure  model 

A  stress-based  failure  initiation  criterion  is  combined 
with  orthotropic  softening  (ductile  failure)  to  find  the 
material  failure  model.  A  complete  definition  of  this 
model  entails  the  knowledge  of  six  (three  normal  and 
three  shear)  failure  initiation  stresses  and  six  corre¬ 
sponding  fracture  energies.  However,  due  to  the  planar- 
isotropic  nature  of  the  material  at  hand,  only  four 
(unique)  failure  initiation  stresses  and  four  fracture 
energies  were  needed.  The  values  for  these  quantities 
were  determined  previously  in  the  multi-length  scale 
computational  analysis  of  the  material  at  hand  [10]. 

3.2.4  Erosion 

To  prevent  run-time  errors  from  occurring  as  a  result 
of  highly  distorted  Lagrangian  cells,  an  erosion  model 
is  also  defined  for  the  material  at  hand.  To  ensure  that 
a  constant  erosion  criterion  is  defined  for  all  the 
materials  used  in  the  present  work,  erosion  is  set  to 
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occur  when  the  instantaneous  geometric  (elastic  + 
plastic  +  crack)  strain  reaches  a  value  of  0.65. 

It  should  be  noted  that  some  of  the  parameters  (e.g. 
A2,  A3,  B0,  Bt  )  were  left  independent  of  the  orientation 
of  MWCNT  mat  orientation  (forest  versus  ply).  This 
was  done  because  the  terms  in  question  had  a  sec¬ 
ondary  effect  on  the  computational  results,  as  con¬ 
firmed  in  the  previous  work  [10].  Since  B0  and  Bi 
were  set  to  0.0,  the  values  of  Cv  and  Trei  are  hence 
immaterial.  Conversely,  model  parameters,  which 
have  a  first-order  effect  (e.g.  failure  strengths  and 
fracture  energies,  Table  3)  were  assessed  and  they 
are  significantly  different  for  the  forest  and  in-plane 
MWCNT  mats. 


3.3  Armour  hard-facing  with  MWCNT-reinforced 
PVEE-matrix  composite  mats 

In  this  section,  the  potential  benefits  of  hard  facing 
the  E-glass  fibre-mat  reinforced  armour  with 
MWCNT-reinforced  PVEE-matrix  composite  mats  is 
explored.  Toward  that  end,  a  monolithic  E-glass/ 
PVEE-matrix  armour  with  a  thickness  of  12.9  mm  is 
selected  as  the  ‘control’.  The  PVEE-matrix  in  this 
armour  is  doped  with  a  0.5-vol%  of  MWCNTs,  follow¬ 
ing  the  recommendation  made  in  the  prior  work  [10] 
with  regard  to  the  benefits  of  doping  the  PVEE  matrix 
with  the  MWCNTs.  The  ballistic-protection  perform¬ 
ance  of  the  control  with  respect  to  a  30  caliber  FSP 
has  been  assessed  by  impacting  the  armour  with 
the  FSP.  The  FSP’s  initial  velocity  is  varied  in  a 
range  between  400  and  425  m/s. 

The  results  of  this  calculation  are  displayed  in 
Figs  3(a)  and  (b).  The  results  displayed  in  Fig.  3(a) 
suggest  that  the  initial  velocity  of  the  projectile, 
which  ensures  armour  penetration  with  a  minimal 
residual  velocity  is  within  the  range  of  405-415  m/s. 
It  should  be  noted  that  the  results  displayed  in 
Fig.  3(a)  and  its  subsequent  counterparts  are  given 
for  a  minimal  time  period  required  to  attain  a  con¬ 
stant  residual  velocity  of  the  projectile  in  each  case. 
While  all  of  the  calculations  carried  out  in  the  present 
work  are  of  a  deterministic  nature,  this  velocity  will 
be  referred  to  in  the  remainder  of  the  manuscript  as 
V50  velocity  (the  initial  velocity  at  which  the  prob¬ 
ability  for  armour  penetration  by  the  projectile  is 
~50  per  cent).  Note  that  a  negative  value  of  the  term¬ 
inal  velocity  indicates  that  the  FSP  has  been  fully 
defeated  by  the  armour  and  it  is  forced  to  bounce 
back.  A  material  deformation/failure  plot  for  the  con¬ 
trol  after  impact  with  an  FSP  (initially  propelled  at  a 
velocity  of  415  m/s)  is  displayed  in  Fig.  3(b).  Regions 
associated  with  different  dominant  deformation/ fail¬ 
ure  modes  are  indicated.  The  results  displayed  in 
Fig.  3(b),  as  well  as  the  results  displayed  in 
Figs  4(b),  5(b),  6(b),  and  7(b)  show  no  evidence  of 


Fig.  3  (a)  Temporal  evolution  of  a  0.30  caliber  FSP 

velocity  during  impact  with  a  12.9-mm  thick 
monolithic  E-glass  continuous  fibre-reinforced 
MWCNT-doped  poly- vinyl-ester-epoxy  matrix 
composite  armour  and  (b)  deformation/ 
damage  distribution  for  the  case  of  the  FSP’s 
initial  velocity  of  415  m/s 

FSP  erosion  despite  the  fact  that  an  erosion  strain 
of  0.65  for  this  material  was  used  in  the  calculations. 
This  value  for  the  FSP  was  obtained  in  the  previous 
work  [13]  in  which  the  interactions  between  an  FSP 
and  ceramic-hard-faced  armour  was  analysed. 

Next,  the  effect  of  hard-facing  the  monolithic 
armour,  discussed  above,  with  a  0.2-mm  thick  PVEE- 
matrix  mat  reinforced  with  10  or  30  vol%  MWCNTs 
aligned  in  a  direction  parallel  to  the  armour  faces 
are  given.  The  results  of  the  ballistic  performance 
analysis  of  the  resulting  hybrid  armour  are  shown  in 
Figs  4(a)  and  (b)  and  5(a)  and  (b)  for  the  cases  of 
10  and  30  vol%  MWCNTs,  respectively.  The  results 
displayed  in  these  figures  indicate  that  the  presence 
of  the  0.2-mm  thick  PVEE-matrix  mat  reinforced 
with  either  10  or  30  vol%  MWCNTs  aligned  in  a  direc¬ 
tion  parallel  to  the  armour  faces  does  not  noticeably 
affect  either  the  V50  (remains  within  the  range  of 
405-415  m/s)  or  the  distribution  of  the  impact 
induced  damage  in  the  armour. 

Finally,  the  effect  of  hard-facing  the  monolithic 
armour,  discussed  above,  with  a  0.2-mm-thick  PVEE- 
matrix  mat  reinforced  with  10  or  30  vol%  MWCNTs 
aligned  orthogonal  to  the  armour  faces  (MWCNT 
forest).  The  results  of  the  ballistic  performance 
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Fig.  4  (a)  Temporal  evolution  of  a  0.30  caliber  FSP 

velocity  during  impact  with  a  12.9-mm  thick 
hybrid  armour  E-glass  continuous  fibre- 
reinforced  MWCNT-doped  PVEE  matrix 
composite  hard-faced  with  a  0.2-mm  thick 
10vol%  MWCNT  ply  mat  and  (b) 
deformation/damage  distribution  for  the  case 
of  the  FSP’s  initial  velocity  of  415  m/s 

analysis  of  the  resulting  hybrid  armour  are  shown  in 
Figs  6(a)  and  (b)  and  7(a)  and  (b)  for  the  cases  of  10 
and  30  vol%  MWCNTs,  respectively.  The  results  dis¬ 
played  in  these  figures  indicate  that  the  presence  of 
the  0.2-mm  thick  PVEE-matrix  mat  reinforced  with 
10  vol%  MWCNT  forests  does  not  noticeably  affect 
either  the  V50  (remains  within  the  range  of  405- 
415  m/s)  or  the  distribution  of  the  impact  induced 
damage  in  the  armour.  Conversely,  in  the  case  of 
30  vol%  MWCNT  forest  mats,  a  noticeable  but 
minor  (5-10  m/s)  increase  in  the  V50  is  observed, 
Fig.  7(a),  whereas  no  significant  changes  to  the  distri¬ 
bution  of  impact  induced  damage  in  the  armour  can 
be  detected,  Fig.  7(b).  It  should  be  noted  that  the 
results  shown  in  Fig.  7(b)  correspond  to  a  post 
impact  time  that  which  is  a  few  per  cent  longer 
than  that  corresponding  to  Fig.  5(b),  which  accounts 
for  the  fact  that  the  extent  of  damage  in  the  armour  is 
somewhat  higher  in  Fig.  7  (b) . 

The  observed  improvement  in  the  ballistic  per¬ 
formance  of  the  E-glass/doped-PVEE-matrix  compo¬ 
site  armour  hard-faced  with  a  layer  of  MWCNT -forest 
mats  is  consistent  with  the  observation  made  by  Cao 
et  al.  [14]  that,  in  a  MWCNT-forest  mat,  it  is  difficult 


Fig.  5  (a)  Temporal  evolution  of  a  0.30  caliber  FSP 

velocity  during  impact  with  a  12.9-mm  thick 
hybrid  armour  E-glass  continuous  fibre- 
reinforced  MWCNT-doped  PVEE  matrix 
composite  hard-faced  with  a  0.2-mm  thick 
30  vol%  MWCNT  ply  mat  and  (b) 
deformation /damage  distribution  for  the  case 
of  the  FSP’s  initial  velocity  of  415  m/s 

for  nanotubes  to  buckle  independently  because  of  the 
proximity  of  the  neighboring  tubes.  The  cooperative 
nature  of  the  buckling  results  in  a  self- organized, 
zigzag-folded  morphology,  which  is  the  most  space- 
efficient  and  energetically  favourable  configuration 
for  huge  numbers  of  nanotubes  to  adopt  under 
large  compressive  strains  brought  about  by  the  FSP. 
The  folding  associated  with  this  zigzag  buckling 
does  not  require  any  extra  space  to  accommodate 
the  vertical  deformations  and  acts  as  a  very  efficient 
mechanism  for  absorption  of  the  kinetic  energy  car¬ 
ried  out  by  the  FSP.  In  addition,  the  upright  orien¬ 
tation  of  the  nanotubes  allows  them  to  make  direct 
contact  with  the  FSP  promoting  its  erosion.  Unfortu¬ 
nately,  the  extent  of  ballistic  protection  improvement 
brought  about  by  the  use  of  30  vol%  MWCNT  forest 
mats  as  armour  hard-facing  observed  in  the  present 
work  is  only  about  1-2  per  cent.  This  finding  is  not 
completely  unexpected  considering  the  fact  that:  (a) 
the  hard-facing  represents  roughly  only  l/65th  of 
the  armour  thickness  and  (b)  due  to  high  compliance 
of  the  armour  in  the  through-the-thickness  direction 
and  low  in- plane  strength  of  the  MWCNT  mats,  the 
FSP  impact  was  found  to  cause  compressive 
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Fig.  6  (a)  Temporal  evolution  of  a  0.30  caliber  FSP 
velocity  during  impact  with  a  12.9-mm  thick 
hybrid  armour  E-glass  continuous  fibre- 
reinforced  MWCNT-doped  PVEE  matrix 
composite  hard-faced  with  a  0.2-mm  thick 
10vol%  MWCNT  forest  mat  and  (b) 
deformation/damage  distribution  for  the  case 
of  the  FSP’s  initial  velocity  of  415  m/s 

deformation  predominately  in  the  E-  glass  PVEE 
matrix  composite.  This  allowed  a  circular  hole  of 
the  size  of  the  FSP  to  be  punched- out  in  the 
MWCNT  forest  mat  without  considerable  mat  com¬ 
pression.  In  other  words,  due  to  the  compliant 
nature  of  the  bulk  section  of  the  armour,  the  full 
potential  of  the  MWCNT  forest  mats  was  not  utilized. 
The  shortcomings  of  the  current  architecture  of  the 
E-glass  PVEE  matrix  composite  armour  hard-faced 
with  MWCNT  forest  mats  can  be  overcome  by:  (a) 
the  use  of  multiple  layers  of  these  MWCNT  forest 
mats  throughout  the  full  thickness  of  the  bulk  section 
of  the  armour;  (b)  depositing  MWCNT  forests  onto 
ceramic  tiles  and  using  such  tiles  for  hard- facing  as 
in  the  case  of  the  multi-functional  composite 
armour  depicted  in  Fig.  1.  In  addition,  following  the 
observations  of  Veedu  et  al.  [11]  that  interlacing  of 
fibre-mat/polymer  matrix  composites  with  MWCNT 
forest  mats  increases  tensile  transverse  properties, 
the  delamination  resistance  of  the  armour  can  be 
perhaps  improved  by  placing  dual-sided  MWCNT 
forest  mats  closer  to  the  back/ spall  face  of  the 
armour  where  significant  delamination  typically 
occurs. 


Fig.  7  (a)  Temporal  evolution  of  a  0.30  caliber  FSP 

velocity  during  impact  with  a  12.9-mm  thick 
hybrid  armour  E-glass  continuous  fibre- 
reinforced  MWCNT-doped  PVEE  matrix 
composite  hard-faced  with  a  0.2-mm  thick 
30  vol%  MWCNT  forest  mat  and  (b) 
deformation/ damage  distribution  for  the  case 
of  the  FSP’s  initial  velocity  of  415  m/s 

The  results  obtained  in  the  present  work  reveal  that 
despite  the  fact  that  all  the  calculations  were  of  a 
deterministic  nature,  generally  no  single  valued  V50 
is  found.  This  is  clearly  seen  in  Fig.  8  in  which  the 
terminal  FSP’s  velocity  is  plotted  as  a  function  of 
the  projectiles  initial  velocity.  The  results  displayed 
in  this  figure  clearly  demonstrate  that  there  is  a 
range  of  initial  FSP  velocity  in  which  the  projectile 
may  or  may  not  penetrate  the  armour  depending  on 
the  exact  value  of  its  initial  velocity.  This  Ending 
demonstrates  the  full  complexity  of  the  projectile/ 
armour  interaction  accompanied  by  plastic  defor¬ 
mation,  failure,  and  erosion  of  the  participating 
materials.  The  oscillations  in  the  terminal  velocity 
observed  in  Fig.  8  cannot  be  readily  explained.  Never¬ 
theless,  the  most  logical  explanation  is  that  the  inter¬ 
actions  between  the  freed  nodes  (produced  after 
erosion  of  the  Lagrangian  elements)  leads  to  a  com¬ 
plex  multi-body  system,  which  can  result  in  different 
terminal  velocity  of  the  FSP. 

Similar  oscillatory  variations  of  the  FSP  projectile 
velocity  with  a  change  in  the  initial  FSP  velocity 
have  been  previously  observed  in  ceramic-hard-faced 
light  composite  armour  [13].  Hence,  this  observation 
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Fig.  8  Variation  of  the  projectile’s  terminal  velocity  as  a 
function  of  its  initial  velocity.  Negative  values  of 
the  terminal  velocity  denote  the  defeat  of  the 
project  by  the  armour 

cannot  be  readily  characterized  as  a  deficiency  of  the 
present  material  model  or  the  computational  pro¬ 
cedure.  Rather,  this  observation  reflects  the  complex 
nature  of  the  FSP/ armour  interactions  where  the 
extent  and  the  distribution  of  damage  in  armour 
regions  surrounding  the  FSP,  and  not  only  the 
region  below  the  projectile,  affect  the  outcome  of  an 
FSP/armour  interaction.  Furthermore,  similar  uncer¬ 
tainties  in  determining  V50  experimentally  are  gener¬ 
ally  encountered. 

Considering  the  fact  that  the  uncertainty  in  the 
determination  of  V50  is  comparable  with  the  highest 
observed  increase  in  this  quantity  with  changes  in 
armour  architecture,  one  must  be  careful  about 
claiming  that  the  use  of  MWCNT  forest  strike  face 
can  indeed  lead  to  improvements  in  the  ballistic  per¬ 
formance  of  light  composite  armour.  This  point  has 
been  partly  confirmed  in  the  ongoing  experimental 
investigations  [15]  of  the  E-glass  mat  reinforced 
PVEE-matrix  armour  hard-faced  with  ca.  17  per  cent 
in-plane  MWCNT-reinforced  PVEE-matrix  compo¬ 
site.  Preliminary  results  showed  no  statistically  sig¬ 
nificant  improvement  in  V50  relative  to  the  same 
armour  without  hard-facing. 


4  SUMMARY  AND  CONCLUSIONS 

Based  on  the  results  obtained  in  the  present  work,  the 
following  summary  and  main  conclusions  can  be 
made. 

1.  A  simple  micro-mechanics  model  for  the  material 
elastic  properties  and  a  lower-bound  estimate  for 


the  strength  and  failure  properties  are  used  to  con¬ 
struct  a  ballistic  material  model  for  MWCNT  forest 
mats  (mats  in  which  MWCNTs  are  aligned  orthog¬ 
onal  to  the  mat  faces)  containing  a 
poly-vinyl-ester-epoxy  matrix. 

2.  A  series  of  transient  non-linear  dynamics  analyses 
of  the  impact  of  E-glass/PVEE-matrix  composite 
armour  hard-faced  with  a  MWCNT-forest  mat 
with  a  fragment  simulating  projectile  is  carried 
out.  The  results  obtained  showed  that  at  10  vol% 
of  MWCNTs  in  the  hard- face  mat,  no  measurable 
improvement  is  obtained.  At  30  vol%  of 
MWCNTs,  a  1-2  per  cent  gain  in  the  ballistic  pro¬ 
tection  performance  of  the  armour  is  observed. 
Although  this  enhancement  is  somewhat  ques¬ 
tionable  considering  the  fact  that  significant 
uncertainty  exists  relative  to  the  determination  of 
V50,  it  appears  that  the  observed  lack  in  ballistic 
resistance  improvement  is  related  to  inferior 
in-plane  strength  properties  of  the  composite 
laminate  armour. 

3.  Since,  similar  improvements  in  the  armour  pen¬ 
etration  resistance  were  not  obtained  in  the  case 
of  armour  hard-facing  with  the  mats  containing 
30  vol%  of  MWCNTs,  but  with  the  MWCNTs 
lying  parallel  with  the  armour  faces,  the  present 
work  revealed  that  the  mode  of  nanotube  defor¬ 
mation  plays  an  important  role  in  the  effectiveness 
of  MWCNT-based  mats  as  armour  hard-facing. 

4.  Potential  directions  for  obtaining  significant 
improvements  in  the  ballistic  protection  perform¬ 
ance  of  the  armour  through  the  use  of  MWCNT 
mats  are  identified  and  discussed. 

5.  The  computational  approach  presented  in  the  pre¬ 
sent  work  showed  that  transient  non-linear 
dynamics  simulations  of  the  projectile/armour 
interactions  can  provide  critical  insight  into  that 
deformation  and  fracture  mechanisms  controlling 
the  ballistic  performance  of  armour.  In  addition, 
the  computational  procedure  reveals  stochastic 
effects  of  the  projectile/armour  interactions 
which  are  commonly  observed  experimentally, 
and  usually  attributed  to  variability  in  armour 
architecture  and  material  properties. 
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